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LONG-TERM  GOALS 

To  develop  and  test  with  laboratory  and  field  data  a  robust  modeling  framework  that  predicts 
hydrodynamics  and  the  fate  of  terrestrial  and  marine  sediment  in  the  heterogeneous  environment. 

OBJECTIVES 

■  Develop  phase/depth-resolving  numerical  models  for  bottom  sediment  transport  and  surf  zone 
hydrodynamics.  Validate  these  models  using  data  measured  in  CROSSTEX  and  other  field 
experiments. 

■  Identify  the  role  of  undertow  current,  steep  waves  and  breaking  wave  turbulence  on  sediment 
suspension  events. 

■  Develop  simplified  phase-resolving  formulations  for  concentrated  sediment  transport,  suspended 
load  transport  and  its  near-bed  boundary  conditions  under  breaking  waves. 

APPROACH 

A  two-phase  flow  model  (Hsu  et  al  2004;  Amoudry  et  al.  2007)  is  extended  to  model  sand  transport 
driven  by  measured  random  wave-current  forcing  during  CROSSTEX  (Scott  et  al.  2007).  Measured 
near-bed  time  series  of  flow  velocities  are  analyzed  to  extract  (approximately)  the  time  series  of 
averaged  wave  velocities  and  turbulence  quantities  (TKE  and  dissipation  rate).  These  time  series  are 
then  used  to  drive  the  two-phase  flow  model  and  to  calculate  and  resulting  sand  transport  from 
immobile  sand  bed  to  the  dilute  region  of  about  5  cm  above  the  immobile  bed.  Physical  experiments 
also  provide  measured  suspended  sand  concentration  in  the  intermediate  to  dilute  region  (volume 
concentration  ~<15%),  which  can  be  used  to  identify  the  suspension  event  and  to  validate  the 
numerical  model.  On  the  other  hand,  two-phase  model  results  for  concentrated  sediment  dynamics  can 
bridge  the  missing  information  at  the  regime  of  concentration  from  about  15%  to  random-close 
packing  (-63%),  which  is  difficult  to  measure  in  the  physical  experiment. 

Because  it  is  non-trivial  to  separate  turbulence  from  the  wave  motion  in  the  measured  velocity  time 
series,  also  because  it  is  perhaps  unclear  at  this  point  what  is  the  major  cause  (e.g.,  steep  wave  or 
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breaking  wave  turbulenee)  of  the  strong  sediment  suspension  event,  it  is  eritieal  to  utilize  a  more 
powerful  data  analysis  tool  than  eommonly  used  Fourier  analysis.  Wavelet  analysis  is  able  to  identify 
not  only  the  steep  oseillation  events  but  also  their  (temporal)  seale.  Henee,  more  information  ean  be 
extraeted  from  the  measured  time-series  when  using  wavelet  analysis  in  eonjunetion  with  Fourier 
analysis.  Wavelet  analysis  on  the  CROSSTEX  data  is  investigated  in  eollaboration  with  Dr.  Nieolas 
Seott  of  Naval  Researeh  Laboratory. 

WORK  COMPLETED 

CROSSTEX  is  a  eollaborative  researeh  projeet  aiming  at  quantifying  various  meehanisms  of  wave- 
and  eurrent-indueed  sediment  transport  in  the  surf  zone.  The  physieal  experiment  of  CROSSTEX  was 
eompleted  in  Sep  2005  at  O.  H.  Hinsdale  Wave  Laboratory  of  Oregon  State  University  (Seott  2006; 
Seott  et  al.  2007).  The  experimental  design  of  CROSSTEX  is  to  mimie  the  field  eondition  of 
onshore/offshore  sandbar  migration  events  observed  during  Duek  94  (Gallagher  et  al.  1998)  so  that 
possible  meehanisms  eontrolling  the  onshore/offshore  sediment  flux  ean  be  earefully  studied. 

Even  in  a  well-eontrolled  laboratory  faeility,  detailed  sediment  transport  in  the  eoneentrated  regime 
near  the  bed  eannot  be  measured  with  high  temporal  and  spatial  resolution.  In  this  study,  we  utilize  a 
two-phase  model,  whieh  resolves  both  the  bedload  and  suspended  load,  to  link  the  missing  information 
in  the  eoneentrated  regime  of  transport.  The  two-phase  model  for  sand  transport  has  been  validated  in 
the  past  two  years  with  published  laboratory  U-tube  results  (Amoudry  et  al.  2007).  This  year,  we  also 

extend  the  model  with  nonlinear  boundary  layer  terms  (i.e.,  approximate  ~  ^  in  the  2D 

two-phase  equations  with  c  the  wave  eelerity)  to  mimie  sediment  transport  under  progressive  waves. 
The  model  is  validated  with  SISTEX99  data  (Dohmen- Janssen  and  Hanes  2004)  in  order  to  investigate 
the  effeet  of  nonlinear  streaming  on  wave-indueed  transport  (Trowbridge  and  Young  1989;  Henderson 
et  al.  2004)  in  eomparison  with  wave  skewness/asymmetry  and  undertow  eurrent  (e.g.,  Gallagher  et  al. 
1998,  Hsu  and  Hanes  2004).  In  the  meantime,  the  two-phase  model  is  also  modified  to  be  driven  by 
measured  wave  veloeity  and  turbulenee  (both  TKE  and  turbulent  dissipation  rate  need  to  be  used)  and 
reproduee  sediment  eoneentration  similar  to  observation  (see  Eig  1). 

A  eollaborative  effort  between  UE/WHOI  and  Dr.  Nieolas  Seott  (Naval  Researeh  Laboratory)  is 
initiated  earlier  this  year  to  analyze  CROSSTEX  data  using  wavelet  transform.  Morlet  wavelet 
transform  has  been  previously  shown  to  be  effeetive  in  identifying  steep  wave  in  open  oeean  (Seott  et 
al.  2005).  Similar  method  is  applied  here  to  the  measured  veloeity  time-series  (50  Hz  sampling  rate). 
Aeeording  to  analysis  done  in  the  prior  years,  the  eross-shore  veloeity  eomponent  has  both  wave  and 
turbulenee  information  in  the  frequeney  range  1-2  Hz,  while  the  along-tank  veloeity  eomponent  has 
elean  turbulenee  information.  This  suggests  that  the  wavelet  analysis  must  be  applied  eurrently  to 
eross-shore  and  along-tank  eomponents  to  obtain  a  better  understanding  of  steep  wave  and  breaking 
wave  turbulenee.  When  plotting  the  wavelet  transform  results  for  veloeities  along  with  measured 
sediment  eoneentration  time  series  at  different  level  above  the  bed,  we  ean  identify  suspension  events 
driven  by  (1)  steep  wave  only  (2)  steep  wave  and  breaking  wave  turbulenee  and  (3)  breaking  wave 
turbulenee  only  (if  oeeurred).  It  is  our  hope  to  eome  up  with  useful  statisties  relating  sand  suspension 
events  with  steep  wave  and  breaking  wave  turbulenee.  This  information  will  be  used  to  drive  the  two- 
phase  model  in  order  to  further  relate  the  hydrodynamie  foreing  to  net  sand  transport  flux. 

In  the  previous  years,  the  2D  RANS  wave  model  (COBRAS,  Lin  and  Liu  1998)  has  been  extended 
with  2"‘^-order  implieit  seheme  to  better  simulate  field/prototype  seale  surf  zone  proeesses.  However,  it 
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still  requires  lots  of  CPU-time  to  conduct  field/prototype  simulation.  In  the  following  year,  we  will 
collaborate  with  Dr.  Losada’s  group  at  Universidad  de  Cantabria,  Spain.  They  have  been  upgrading 
COBRAS  using  Fortran  90  with  C++  library  and  the  resulting  code  appears  to  compute  several  times 
faster  than  the  our  current  version  of  COBRAS  provided  by  Dr.  Liu  of  Cornell  University. 

RESULTS 

Measured  near-bed  flow  velocity  time  series  during  CROSSTEX  are  used  to  drive  the  two-phase 
model.  The  model  then  reproduces  fluid  velocity,  sediment  concentration,  velocity  and  sediment  flux 
profiles  from  the  immobile  porous  bed  to  the  dilute  regime  of  transport.  For  a  given  wave-current 
forcing,  the  breaking  wave  turbulence  appears  to  enhance  the  magnitude  of  transport.  Based  on  limited 
model  simulations  conducted  at  this  point,  such  transport  is  often  offshore  directed  due  to  strong 
undertow  current:  Concurrent  measurements  of  wave  velocities,  turbulent  velocity  fluctuations  and 
sand  concentration  are  collected  under  random  breaking  wave  trains  near  the  bar  crest  during 
CROSSTEX  (Scott  et  al.  2007).  According  to  measured  data,  strong  sand  suspension  events  are  often 
observed  when  breaking  wave  turbulence  penetrates  into  the  water  column  and  approaches  the  wave 
boundary  layer  (WBL).  Conventional  boundary  layer  sediment  transport  models  cannot  incorporate  the 
effects  of  free-surface  generated  turbulence  on  WBL.  Here,  we  propose  a  new  method  where  both  the 
measured  time-series  of  wave-current  velocity  and  turbulence  quantities  (turbulent  kinetic  energy 
(TKE)  and  dissipation  rate)  are  concurrently  used  to  drive  the  two-phase  model.  Fig  1  presents  the 
model  results  driven  by  measured  80  sec  time  series.  The  lower  three  panels  show  the  model  results  for 
velocity,  turbulence  intensity  and  sand  concentration  profiles  at  three  instants  during  the  passage  of  a 
breaking  wave  train  (top  panel).  When  breaking  wave  turbulence  does  not  approach  the  WBL  (tl), 
model  results  obtained  by  using  measured  turbulence  as  top  forcing  (solid  curve)  and  using 
conventional  approach  (dashed)  are  almost  identical.  Both  methods  predict  similar  sand  concentration 
compared  to  measured  data  (crosses).  However,  when  the  breaking-wave  turbulence  is  strong  at  (t2) 
and  (t3),  model  results  obtained  by  also  using  measure  turbulence  as  forcing  at  top  boundary  shows 
strong  TKE  across  the  WBL.  In  addition,  model  predicted  sand  concentrations  (solid)  are  similar  to 
measured  data  (crosses)  while  conventional  method  fails  to  predict  the  suspension  event  (dashed). 
Despite  the  individual  suspension  events  are  strongly  affected  by  the  breaking  wave  turbulence,  the 
time-averaged  net  transport  rate  for  this  case  is  only  enhanced  by  30%.  This  is  possibly  because  the 
suspension  events  are  usually  of  short  duration  compared  to  the  overall  length  of  the  wave  train. 

On  the  other  hand,  using  time  series  measured  at  other  times,  the  net  transport  rate  can  be  enhanced  by 
factor  2  (not  shown  here).  We  also  find  that  the  result  of  net  transport  rate  is  rather  sensitive  to  how  we 
extract  turbulence  from  the  original  velocity  signal.  Therefore,  we  decide  to  focus  on  more 
sophisticated  data  analysis  technique  using  wavelet  transform  in  order  to  identify  the  suspension  events 
and  their  responses  to  steep  wave  and  breaking  wave  turbulence.  Such  analysis  will  allows  us  to  more 
objectively  quantify  the  occurrence  frequency  of  these  suspension  events  and  their  relation  to  the 
forcing.  Results  obtain  from  such  analysis  can  provide  a  more  clear  guideline  for  numerical  modeling 
of  sediment  transport  under  breaking  waves. 

IMPACT/APPLICATIONS 

The  present  research  efforts  focus  on  developing  and  validating  detailed  numerical  models  for 
sediment  transport  and  wave  hydrodynamics.  We  have  focused  on  refining  the  numerical  schemes  in 
the  two-phase  model  and  RANS  wave  model  so  that  they  will  be  more  robust  and  user-friendly  in  the 
future.  It  expected  that  by  the  end  of  this  4-year  project,  the  two-phase  model  will  be  made  available  to 
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the  research  community.  Recent  effort  using  wavelet  analysis  to  study  sediment  transport  under 
breaking  wave  is  also  new  to  the  coastal  sediment  transport  community. 


RELATED  PROJECTS 


Hsu  is  currently  supported  by  a  NSF-OCE  project:  3D  Multiphase  Sediment  Transport  Modeling 
Framework.  This  project  focuses  on  extending  the  existing  IDV  ensemble-averaged  two-phase  model 
for  more  detailed  3D  eddy-resolving  multiphase  model.  This  model  development  effort  will  be 
extremely  useful  for  the  present  research  efforts  to  understand  and  parameterize  sediment  transport 
under  breaking  waves. 
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Fig  1:  Two-phase  model  is  driven  by  measured  near-bed  wave-current  velocity  and  turbulence 
at  5cm  above  the  bed  during  CROSSTEX  (upper  panel).  In  the  lower  panels  (tl),  (t2)  and  (t3) 
are  model  results  for  fluid  velocity,  turbulent  intensity  and  sediment  concentration.  Black- 
dashed  curves  are  results  without  considering  the  breaking  wave  turbulence  while  the  red 
curves  results  when  measured  breaking  wave  turbulence  is  considered.  Cross  symbols  are 
measured  concentration.  When  the  breaking  wave  turbulence  is  considered,  the  model  is  able 
to  predict  the  suspension  events  observed  during  the  experiment,  suggesting  the  importance  of 
incorporating  breaking  wave  turbulence  in  modeling  the 
bottom  boundary  layer  sediment  transport. 
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